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Abstract 

Many models of supersymmetry breaking, in the context of either supergrav- 
ity or superstring theories, predict the presence of particles with weak scale 
masses and Planck-suppressed couplings. Typical examples are the scalar 
moduli and the gravitino. Excessive production of such particles in the early 
Universe destroys the successful predictions of nucleosynthesis. In particu- 
lar, the thermal production of these relics after inflation leads to a bound on 
the reheating temperature, Tjih ~ 10^ GeV. In this paper we show that the 
non-thermal generation of these dangerous relics may be much more efficient 
than the thermal production after inflation. Scalar moduli fields may be co- 
piously created by the classical gravitational effects on the vacuum state. 
Consequently, the new upper bound on the reheating temperature is shown 
to be, in some cases, as low as 100 GeV. We also study the non-thermal 
production of gravitinos in the early Universe, which can be extremely effi- 
cient and overcome the thermal production by several orders of magnitude, 
in realistic supersymmetric inflationary models. 

July 1999 



*E-mail: Gian.GiudiceOcern. ch 
^E-mail: Igor .TkachevScern. ch 
■tEmail: r iotto®nxth04 . cern . ch 



1 Introduction 



One of the problems facing Planck-scale physics is the lack of predictivity for low-energy 
phenomena. Nonetheless, some information on high-energy physics may be inferred 
indirectly through its effects on the cosmological evolution of the early Universe. If we 
require that the successful predictions of big-bang nucleosynthesis are not significantly 
modified and that the energy density of stable particles does not overdose the Universe, 
severe constraints are imposed on the properties of a large class of supersymmetric 
theories. 

In = 1 supergravity models supersymmetry is broken in some hidden sector 
and gravitational-strength interactions communicate the breaking down to the visible 
sector. In these models there often exist scalar and fermionic fields with masses of the 
order of the weak scale and gravitational-strength couplings to ordinary matter. In 
the following, we will generically refer to them as gravitational relics X. If produced 
in the early Universe, such quanta will behave like nonrelativistic matter and decay at 
very late times, eventually dominating the energy of the Universe until it is too late 
for nucleosynthesis to occur (in the case of scalar fields, long-lived coherent oscillations 
with large amplitudes of the zero mode can pose the same problem) . Typical examples 
of gravitational relics are the spin-3/2 gravitino - the supersymmetric partner of the 
graviton - and the moduli fields - the quanta of the scalar fields which parametrize 
supersymmetric fiat directions in moduli space and seem almost ubiquitous in string 
theory. 

In string models massless fields exist in all known ground states and parametrize 
the continuous vacuum degeneracies characteristic of supersymmetric theories. Possible 
examples of gravitational relics in string theory are the dilaton whose vacuum expecta- 
tion value parametrizes the strength of the gauge forces, and the massless scalar (and 
fermionic superpartner) gauge singlets parametrizing the size of the compactified di- 
mensions. In general, these fields are massless to all orders in perturbation theory in 
the exact supersymmetric limit and become massive either because of non-perturbative 
effects or because of supersymmetry-breaking contributions. In the usual scenarios in 
which supersymmetry breaking occurs at an intermediate scale Ms = 0{^/MpMw), 
moduli acquire masses of the order of the weak scale and, because of their long lifetime, 
pose a serious cosmological problem IQ. The problem becomes even more severe in 
models in which supersymmetry is broken at lower scales, since the moduli are lighter 
than Mw and their lifetime becomes exceedingly long. On the other hand, if supersym- 
metry breaking occurs at a scale larger than Ms and its effect on the observable sector 
is screened, the cosmological problem of gravitational relics is relaxed. In particular. 
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this is the case of the recently-proposed models with anomaly-mediated supersymmetry- 
breaking 0, in which moduli and gravitinos can avoid cosmological difficulties 

Gravitational relics X can be produced in the early Universe because of thermal 
scatterings in the plasma. The slow decay rate of the X-particles is the essential source 
of cosmological problems because the decay products of these relics will destroy the "^He 
and D nuclei by photodissociation, and thus successful nucleosynthesis predictions |§, 
The most stringent bound comes from the resulting overproduction of D + ^He; this 
requires that the X-abundance relative to the entropy density at the time of reheating 
after inflation should satisfy 

— < 10-12. (1) 

s 

(The exact bound depends upon the X mass, and here we have assumed Mx ~ TeV.) 

Neglecting any initial number density, the Boltzmann equation for the number den- 
sity of gravitational relics during the thermalization stage after inflation is 

^ + 3Hnx ^ {(rtot\v\)nl^u, (2) 

where cxtot oc 1/Mp is the total production cross section and riught ~ represents 
the number density of light particles in the thermal bath. The number density nx at 
thermalization is obtained by solving eq. (0), 

— ^ 10-^ (3) 

Comparing eqs. (|l]) and (H), one obtains an upper bound on the reheating temperature 
after inflation 

Tbh < (10^ - 10^) GeV. (4) 

If Trh ~ ^GUT ~ 10^^ GeV, gravitational relics are inconsistent with nucleosynthesis. 
Moreover, if the initial state after inflation is free from gravitational relics, the reheating 
temperature in eq. (H) is too low to allow for the creation of superheavy GUT bosons 
that can eventually produce the baryon asymmetry ||^. 

The goal of this paper is to point out that gravitational relics can be created with 
dangerously large abundances by non-thermal effects. 

One possibility is represented by the non-thermal effect of the classical gravitational 
background on the vacuum state during the evolution of the Universe. The particle cre- 
ation mechanism is similar to the inflationary generation of gravitational perturbations 
that seed the formation of large scale structures. However, the quantum generation 
of energy density fluctuations from inflation is associated with the inflaton fleld which 
dominated the mass density of the Universe, and not a generic, sub-dominant scalar or 
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fermionic field. Gravitational particle creation has recently been employed to generate 
non-thermal populations of very massive particles |T^. In particular, the desired 
abundance of superheavy dark matter particles may be generated during the transition 
from the inflationary phase to a matter/radiation dominated phase as the result of the 
expansion of the background spacetime acting on vacuum quantum fluctuations of the 
dark matter field. Contrary to the generation of density perturbations, this mechanism 
contributes to the homogeneous background energy density that drives the cosmic ex- 
pansion, and is essentially the familiar "particle production" effect of relativistic field 
theory in external fields. 

Another possibility is represented by the non-thermal effects occuring right after 
inflation because of the rapid oscillations of the infiaton field(s). Gravitino production 
is an interesting example of this kind. 

The paper is organized as follows. In sect. 2 we show that the gravitational pro- 
duction of scalar moduli which are minimally coupled to gravity may be much more 
efficient than the one activated by thermal scatterings during reheating. Our findings 
lead to extremely stringent constraints on the reheating temperature after inflation. In 
sect. 3 we study the non-thermal production of gravitinos in the early Universe, showing 
that the helicity-1/2 part of the gravitino can be efficiently excited during the evolu- 
tion of the Universe. This leads to a copious and dangerous generation of gravitinos 



in realistic supersymmetric models of inflation. Previous studies [|TT| of non-thermal 
production of gravitinos have considered only the subleading effects of the helicity-3/2 
components. However, a thorough study of the complete gravitino equations has very 
recently appeared . Although the derivation of the gravitino equations presented in 



our paper is different from that of ref. [T^, our final results agree fully with the analysis 
of ref. 1 12]. Finally, sect. 4 contains our conclusions. 



2 Gravitational Production of Scalar Moduli 

Let us describe here the basic physics underlying the mechanism of gravitational scalar 
particle production. 

We start by canonically quantizing the action of a generic scalar massive field X 
which, at the end, will be identified with a modulus field. In the system of coordinates 
where the line element is given by ds^ = dt^ — a?'{t)dx^, the action is 

S = jdtj d\'^ 

where R is the Ricci scalar. After transforming to conformal time coordinate, where 




{VXf 



(5) 
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the line element is ds^ = a?{ri){drf — doF), we use the field mode expansion 

^ = / i2n)VMv) ["^^^^^^'^'■' + aihlir,)e-''-'] ■ (6) 

Since the creation and annihilation operators obey the commutator relations [a^^, alj = 
S^^\ki — k2)i we obtain the normalization condition /ifc/i^ — hkh], = i (henceforth, the 
dots over functions denote derivatives with respect to r]). The resulting mode equation 
is 

hkiv)+^liv)hkiv)=0, (7) 

where 

oul = e + M],a^ + {6^ - 1)- . (8) 

a 

The parameter ^ is 1/6 for conformal coupling and for minimal coupling. The equation 
of massless conformally coupled quanta reduces to the equation in fiat space-time and, 
in this case, there is no particle production. When Mx 7^ or ^ 7^ 1/6, conformal 
invariance is broken and particles are created. In the following we will be especially 
interested in the case of scalar particles with non-conformal coupling. In particular, 
the case of minimal coupling {C, = 0) is certainly of physical interest. Just to give one 
example, in the toroidal compactification of the type IIB string theory, the modulus field 
describing the volume of the internal compactified dimensions is minimally coupled to 
gravity in the four dimensional effective action. 

The differential equation (0) can be solved once the boundary conditions are sup- 
plied. Since the annihilation operator is just a coefficient of an expansion in a particular 
basis, fixing the boundary conditions is equivalent to fixing the vacuum. We choose the 
initial conditions as 

hkiVo) = ^k^''^, hkivo) = -iujhk{r]o), (9) 
corresponding to a vanishing particle density at rj = tjq. To obtain the number density 
of the produced particles, we perform a Bogolyubov transformation from the vacuum 
mode solution with the boundary condition aX rj = rjQ into the one with the boundary 
condition at rj = r]i (any later time at which the particles are no longer being created). 
Defining the Bogolyubov transformation as hP^ij]) = akhPj^irf) + Pk^l^^iv) (the super- 
scripts denote where the boundary condition is set), we obtain the following number 
density of produced particles: 

1 1'°° 

^^(^1) = ^T^^T^ / dkkW. (10) 
2TT^a-^[rii) Jo 

One should note that the number operator is defined at 771 while the quantum state 
(approximated to be the vacuum state) defined at r/o does not change in time in the 
Heisenberg representation. 
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Let us now discuss the structure of the mass term of the X scalar particle. The flat 
directions in the supersymmetric potential, corresponding to massless X particles, are 
lifted by supersymmetry breaking and nonrenormalizable effects in the superpotential. 
An important observation for us is that in the early Universe global supersymmetry 
is broken (for instance because of the false vacuum energy density of the infiaton). In 
supergravity theories, supersymmetry breaking is transmitted by gravitational interac- 
tions and the supersymmetry-breaking mass squared is naturally ChH^, where H is 
the Hubble parameter and Ch = C>{1). To illustrate this effect, consider a term in the 
Kahler potential of the form 



where I is the field which dominates the energy density p of the Universe, that is 
p — {J d'^OPl). During inflation, / is identified with the infiaton field and p = V{I) = 
SH^Mp. The term (|ll]) therefore provides an effective X mass M]^ = 3ChH^- The 
same mass term will appear during the coherent oscillations of the infiaton field / after 
inflation. This example can be easily generalized, and the resulting potential of the 
scalar moduli fields during inflation is of the form V = H'^Mp V(|X|/Mp). We conclude 
that moduli quanta will be gravitationally produced in the early Universe with a mass 
squared of the form 



Here mx accounts for the mass term generated by any possible source of supersymmetry 
breaking whose F-term is not dominating the energy density of the Universe during 
inflation, but persists in the zero-temperature limit. For this reason, we expect mx 
to be of the order of TeV, the present effective supersymmetry-breaking scale. It is 
important to bear in mind - though - that the case Ch = is certainly a possibility. 
Indeed, in supergravity models which possess a Heisenberg symmetry, supersymmetry 
breaking makes no contribution to scalar masses, leaving supersymmetric flat directions 
unmodified at tree-level. No-scale supergravity of the SU{N, 1) type and the untwisted 
sectors from orbifold compactifications are special cases of this general set of models. 

We have numerically integrated eq. (J^), in presence of an infiaton with quadratic 
potential, up to the epoch when the X field starts oscillating on all scales and the 
field fluctuations are transformed into particles which redshift as thereafter. This 
procedure determines a well-defined quantity, the total number of particles nxo,^, which 
is constant at late times and is easy to compute numerically. However, the calculation is 
technically difficult in the regime of very small masses. In this regime a better strategy 
is represented by adopting observables which are independent of mx at small mx- In 
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this limit - and Ch ^ - the present-day energy density of produced X-particles is an 
appropriate quantity, since it is independent of rrix (see ref. [0] for a complete discussion 
of this point). Therefore, for small masses, the number density nx is better inferred from 
the ratio px/^x- On the other hand, we expect that at Ch ~ 1 the number density of 
produced X-particles will be independent of mx, since the field fluctuations will enter 
the oscillating regime at all scales at the epoch in which ^ m^. Therefore, in the 
limit of relatively large Ch-, it is more convenient to use nx- 

In our numerical results, we have normalized nx by the entropy density. In doing 
so, we have assumed that the reheating temperature after inflation is Trh = 10^ GeV 
and the mass of the inflaton field is m/ = 10^^ GeV. The results for scalar fields with 
minimal coupling = 0) are summarized in fig. |l|. This figure shows the limiting 
behaviour of nx/s as discussed above, which allows to extract the relevant physical 
quantities at arbitrarily small mx- In particular, notice that for small Ch-, the ratio 
nx/s scales like ni^ . If the reader wishes to adopt different values for the reheating 
temperature or the inflaton mass, the data in fig. |I| have to be multiplied by the ratio 

( ^ ns) 

109 GeVy U013 GeVy ■ ^ ' 

Furthermore, if one wants to include the effect of significant entropy release at some 
late time after reheating, the data in fig. |l| have to be divided by the amount of entropy 
increase in the comoving volume, 7 = Sfinai/smitiai- 

From fig. we infer that scalar moduli quanta are very efficiently produced by 
gravitational effects. In particular, for the realistic case mx ~ ~ 1 TeV (or 

mx/nii ~ 10~i°), the ratio nx/s turns out to be extremely large for small Ch 

^ ^ 5 X 10-^ (14) 

This result is seven orders of magnitudes above the limit in eq. ([1|). In this case, the non- 
thermal production is so efficient that the thermal scatterings during reheating become 
completely irrelevant and the upper bound on the reheating temperature in order to get 
nx/s ^ 10~i^ becomes as low as 100 GeV. Only when the parameter Ch approaches 
unity, regardless of the value of mx, does the number density of X particles become 
acceptably small. 

Another way of presenting the bound in eq. (|I]) is to compute the energy density in 
scalar moduli particles normalized to the critical density pc in the Universe, fix = Px/Pc- 
Were these particles stable, the parameter Qx would be 

.,,,.4X10"(^)(!^). ,15) 
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Figure 1: The ratio nx/s as a function of Ch for scalar moduli particles X minimally 
coupled to gravity and with squared masses = m\ + ChH"^- In the figure nix is 
expressed in units of the inflaton mass m/, and Trh = 10^ GeV, mj = 10^^ GeV. 

where h is the Hubble rate in units of 100 km sec~^ Mpc~^. Therefore the bound in 
eq. (|l|) can be recast into the form 



nxh'^ < 4 X 10" 



/ mx 
VTeV 



(16) 



The parameter Qx is shown in fig. 2. Again, this figure was calculated assuming 
Trh = 10^ GeV and mj = 10^^ GeV. If one wishes to adopt different values of Tjih 



and rrij or to take into account some entropy release, the data in this figure have to be 
multiplied by the factor 
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RH 



Vl09 GeV; V1013 GeV; 



(17) 



Figures 3 and 4 represent the ratio nx/s and Qxh"^ as a function of ^ for Ch = 0. 
From the data one can infer that scalar moduli particle production is extremely dan- 
gerous unless ^ is very close to the conformal value ^ = 1/6, where particle production 
shuts off for small mx- 

Our findings are also relevant for the idea of enhanced symmetries and the ground 



state of string theory |]13[. We know that there are often points in the moduli space with 
maximally enhanced symmetry, i.e. points where all of the moduli are charged under 
some symmetry (which may be continuous or discrete). These moduli do not suffer from 
the cosmological moduli problem generated by their large coherent oscillations. For 
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Figure 2: Ratio of the energy density in scalar moduli particles minimally coupled to 
gravity to the critical density, as a function of Ch- Conventions are the same as in fig. 
1. 



m 



10-11 
10-12 

10-13 

10- 
10- 
10- 
10- 
10- 
10- 

10-20 



-\ \ 




1 ' : 


\ 


\ \, 
\ \, 
\ \. 
X \, 






\ \^ 






X \ 






X \ 






X \ 






X \ 






X \ 






X \ 






• ^ \ 












\ 






\\ 












v-. /-•-' 






\ 






1 V/ , , 


1 . " 



0.1 



0.2 



0.3 



Figure 3: The ratio nx/s of scalar moduli particles as a function of the parameter ^, for 
Ch = 0. The solid, long-dashed and short-dashed lines correspond to mx/Tni — 0.001, 
mx/mi = 0.01 and mx/mj = 0.1, respectively. 
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Figure 4: Ratio of the energy density in scalar moduli particles to the critical density as 
a function of ^, for Ch = 0. The solid, long-dashed and short-dashed lines correspond 
to mx/frij = 0.001, mx/frii = 0.01 and mx/frii = 0.1, respectively. 



ordinary moduli, finite temperature and/or curvature effects (say during inflation) are 
likely to leave the Universe in a state which does not correspond to the present minimum 
of the moduli potential. The Universe remains in this state until the Hubble parameter 
is of order of the curvature of the potential, after which the system oscillates. The 
moduli typically dominate the energy density of the Universe when they decay, leading 
to catastrophic consequences. However, if the vacuum is a point of maximally enhanced 
symmetry, it is quite natural for the Universe to start out in this state. For example, 
during inflation, even though the potential for the moduli is modified, it can be expected 
that the system remains in the symmetric state. Finite temperature effects also tend 
to prefer states of higher symmetry |jl3| . However, this selection rule does not prevent 
the scalar quanta from being copiously produced by the gravitational effects, unless 
the effective mass of the scalar excitations around the point of enhanced symmetry is 
much larger than the Hubble rate during inflation. Therefore, gravitational production 
of these special moduli fields has to be checked case by case, in order to assess the 
viability of a model. 



3 Non-Thermal Production of Gravitinos 



In this section we compute the non-thermal production of gravitinos in time-dependent 
gravitational backgrounds. Before solving the relevant equations of motion in curved 



space, it is useful to recall the known results of gravitino propagation in flat space. 

The free propagation of massive gravitinos in Minkowski space is described by four 
Majorana spinors ip"", satisfying the Rarita-Schwinger equation 



R" = e^'"'%-fbdcipd + 2m(T''Vf. = 0. (18) 

Here flat space indices are denoted by Latin letters and are contracted by the metric 
1] = diag(+, —,—,—); we also use the convention 

e'''' = +l, a'^' = lba,lb]. (19) 

The mass term in eq. (p!8D explicitly breaks supersymmetry. However, we are implicitly 
assuming that supersymmetry is spontaneously broken and eq. (|18|) follows from choos- 
ing a supersymmetric gauge such that the Goldstino field is fully eliminated from the 
Lagrangian. This choice is analogous to the unitary gauge in Yang-Mills gauge theories. 
Using the identity 

e 757fc = 2V^^^ -777 J, (20) 

eq. ([T8|) becomes 

= ^ (7" ^ 7 ■ ^ - 7' ^ 7" ^b) + m (7" 7 ■ - V^") = 0. (21) 



From eq. (pll) we obtain the following two constraints 



J- R = 2i J ■ ip - d ■ ip) + 3m J ■ ip = 0, (22) 
d-R = m(^7-V-5-V') = 0. (23) 

For m 7^ 0, these constraints imply the conditions 'y-ijj = and d-ip = 0, which eliminate 
two spinors out of ip"', leaving the appropriate degrees of freedom for the propagation of 
the spin ±3/2 and ±1/2 components. Replacing eqs. (P^)-(PB|) in eq. (^), we obtain 
that the propagation of physical states obeys the ordinary Dirac equation 

{i^ -w^tp" = 0. (24) 

For m = 0, eq. (|18D possesses a gauge symmetry 5^ipa = dae, and we can choose a 



gauge fixing such that 7 ■■?/' = 0. Then, eq. (p2[) leads to d ■ ip = 0. However, the choice 



J ■ ip = still allows gauge transformation subject to = 0, and a further condition 
is required to fully fix the gauge. This condition eliminates an additional fermionic 
component, leaving only the ±3/2 spin states as physical particles. Notice the complete 
analogy with the case of spin-one particles, in which the Lorentz condition d ■ A = 
follows from the equation of motion in the massive case and it is a gauge choice in the 
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massless case. Also, for a massless spin-one particle, the gauge fixing requires a further 
condition, which eliminates the scalar degree of freedom. 

Let us turn now to the case of curved space, in which the gravitino equation of 
motion becomes 

^ = 0. (25) 

Greek letters denote space-time indices and Latin letters denote tangent-space indices. 
Gamma matrices and the Levi-Civita symbol with curved indices are defined by 

r ^ e^^ e'^^"'^ ^ ee^^ele^^e^e'^"^', (26) 

where is the vierbein. is the covariant derivative with respect to the spinorial 
structure modified to also reproduce the mass term 

V,, = D^ + ^m%, D^ = d^ + ]^u}^abcr^\ (27) 

The spin connection uo^ab can be obtained (in first-order formalism) by solving the su- 
pergravity equation of motion where the gravitino, the vierbein, and the spin connection 
are treated as independent fields. One finds (see e.g. ref. |T^) 



^^lab — 2 i~CfJ.ab + Cabfi + Cbfia) , (28) 

= d,et-d.e;-^^,r^,, (29) 

where ip^ = 'ipj^'j^. 

For our cosmological considerations, we are interested in the case of spatially fiat 
Friedmann- Robertson- Walker metrics, in which the line element can be written as ds"^ = 
o?{ri){diff — dx^). Here x = (x^,x^,x^) are comoving space coordinates, x° = 77 is the 
conformal time and a is the scale factor, such that = drj/dt. Thus, the vierbein and 
the metric can be written as = a5jj, = a~M^, g^y = a^rj^^, gf^" = a~'^ri^'' . In this 
case, the covariant derivative in eq. (pTj) reduces to 



= + £1 (7m7° - 7°7m) + (30) 

where the dot denotes the derivative with respect to conformal time. Here we have 
dropped from eq. ( pOf ) the torsion term, bilinear in the gravitino field. This term is 
crucial for establishing the consistency of supergravity However, for our consid- 
erations, it can be safely ignored since the number density of produced gravitinos is 
sufficiently small. During the cosmological epochs we will consider, both m and H are 
much smaller than the Planck mass. 
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Using the identity in eq. (|20|), we can rewrite eq. ( pS] ) in the Friedmann- Robertson- 
Walker metric as 



+ {^^^^^ ~ 2^^^° + ^Yi^i ■ ^ - g^^i ■i^)+m . - ^/^) = 0, (31) 

d^ip^, = d^ip^ + 2a/ail)^ . The condition V ■ R = gives the 



where ^ = 7^9^ and d -ip — ly 
following algebraic constraint^: 



2=1 

1 



3a(f| + 



a a 
-2^ + ^ 



- 3m^ 7" + 2z- 



.m 



(33) 
(34) 



Here we have also considered the possibility that the gravitino mass is time dependent, 
since in general it is determined by time-varying supersymmetry-breaking background 
fields. In the flat space-time limit with constant gravitino mass, one gets c = —7° and 
recovers the gauge flxing condition 7 ■■?/' = 0. 

Since spatial translations generate an exact symmetry of space-time, it is convenient 
to expand the gravitino fleld in momentum modes, iIj^{7],x) ~ / d^ke~'^^'^il)i{ri). In the 
following we will consider the equation of motion for a single momentum mode and 
choose the coordinates such that k is along the direction. For simplicity, we drop the 
index k from ip'i. Because of the antisymmetric properties of the Levi-Civita symbol, the 
equation i?*^ = does not contain time derivatives and therefore describes an algebraic 
constraint on the gravitino momentum modes, which is given by 



d 



{d-f) (71^1+72^2), 
k 



(35) 
(36) 



where k = \k\. It is convenient to describe the remaining independent flelds with two 
Majorana spinors il)i/2 and V'3/2; deflned such that 



C 73 d 7i ?/'i/2. 



(37) 



^The parameter c can be expressed in terms of the background energy density p and pressure p 
using the Einstein equations 



c = 



(32) 



This expression agrees with the result found in ref. ui 
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= 71 + ^3/2, 

= ^'^^ (t! - ^ ^3/2^ 

= J| (c?- 7^)71 V'l/2- 



(38) 
(39) 

(40) 



We can show that, in the flat hmit and on mass-shell, ipi/2 and 1^3/2 correspond to the 
±1/2 and ±3/2 helicity states by explicitly constructing the helicity projectors. Since 
the gravitino field is built from the direct product of spin 1/2 (spinorial indices) and 
spin 1 (vectorial indices) states, the helicity projectors can be decomposed using the 
appropriate Clebsch-Gordan coefficients 



^±3/2 - -P±l/2-P±l 
/I 



' ±1/2 



(41) 
(42) 



The helicity projectors acting on spinorial {P±i/2) and vectorial {P±io) indices are 

P±i/2 = i(l±*7V), 



^±1 = 73(0, ±1, 1, 0), Pi; = 0, 0, vw+^) 



(43) 
(44) 



It is now easy to verify that, in the fiat limit and on mass-shell, eqs. ( |1T| ) and ( |1^ ) 
project onto ipi/2 and V'3/2; respectively. One can also verify that the normalization 
chosen in eqs. (^7|)-(^) insures that the Rarita-Schwinger Lagrangian leads to canonical 
kinetic terms for the fields V'1/2 and 'ip^/'^ 

The equations of motion for the fields ipi/2 and ipz/2 are derived from eq. (^), 

^"3/2 = 0, 



i7°9o ± i — 7° — ma + k^y^ 
2a 

z7°9o ± z— 7° -ma + k(^A + iBj^^ 7' 



A 



B 



2a 
1 



3(f + 



2m 



3 f^±m^ 



m 



aa a 



+ 



1pl/2 = 0, 



4m — ± 3m — 4— mm 



a m 



2^ - — ±3m^— ± 



a^ a^ 



m — 



a^ ma 



■^Notice that the Unear combination 



(46) 
(47) 
(48) 

(49) 
(45) 



defines the ±1/2 hehcity state adopted in ref. Q. The states V1/2 and Si7iV'* differ by a time- 
dependent function. 
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With the field redefinition in eqs. (|37|)-(^0D, the massive Rarita-Schwinger Lagrangian 
has been diagonahzed. The ±3/2-hehcity states satisfy the same equation of motion as 
an ordinary Dirac particle, except for a 5/2 coefficient replacing the usual 3/2 in front 
of the a/ a term. This coefficient is determined by the field scaling with a, and it can 
be simply obtained by recalling that the Lagrangian density should scale as an energy 



density, £ ~ a Since dn ~ aP and 



a simple inspection of the kinetic terms 



and ■il)^' 



~ a 



^5/2 for 



in the corresponding Lagrangians shows that cf) ~ a^^, ip a 
spin-0, spin-1/2, and spin-3/2 particles, respectively. 

The ±l/2-helicity states, on the other hand, satisfy the more complicated evolution 
described by eq. (^7]) . In the de Sitter limit (a = 2a^/a) with constant gravitino 
mass (m = 0), the coefficients in eq. (|47|) become A = cos 6 and B = sin 6, with 
tan(^/2) = a/{ma^). In general, however, the time dependendence of the gravitino 
mass m is related to the evolution of the scale factor. Let us consider the simple case 
of a single chiral superfield $ with minimal kinetic terms. The diagonal time and space 
components of the Einstein equation become 



3M2 



dt 



d^ 



dt 



(50) 
(51) 



Using the expression for the gravitino mass m in terms of the superpotential W , 



Ml ' 



we can write the scalar potential V as 



V = e^'^ 



M'ii 



Ml 



?>\w\- 

~Mf 



m^Ml 



mMp 



am 



dt 



-3 



Replacing eqs. ( pO[ ) and (|51| ) in eq. (|53|), we obtain 



9 d ( 

m= 7 + - -r 

a* a \ a* 



- 3m^ + 2— + 6— 



(52) 



(53) 



(54) 



When this expression for rh is used in eqs. (^oj) and (^9|), we obtain 



(55) 



This condition plays an important role in the solution of the equation for the ±1/2- 
helicity states. 
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To proceed, we explicitly write the momentum expansion for ■?/'i/2 in terms of creation 
and annihilation operators 



X] = a 



(27r)3/2 



r=l,2 



(56) 



where v'^{ri^ k) = u^'-^{ri, —k). For simplicity, we now focus on a single polarization state 
and omit the index r in the following. Choosing a gamma-matrix representation such 
that 

■'-'I 7'^f«, (57, 



1 
-1 



7 = 1 n 11 ) 7 

the spinor u^{ri, k) satisfies the equation of motion (^Tf), i.e 

u± = ^imau± + k{iA =i= B)u^, 



{51 



where u'" = The spinor v{ri,k) satisfies the same equation of motion. The 

treatment of is analogous and it is obtained by just setting A = 1 and i? = in 
the previous equations. 

We can now reduce the system (|58|) of first-order differential equations into a second- 
order differential equation for the function / = G^^^'^u^ with G = A + iB: 



f + 



k^ + vr?Ci' -|- i(ma)' — i—ma + 



G 



G 3G^ 



/ = o. 



G "^ ' 2G 4G\ 
Here we have made use of the property \G\ = 1, see eq. (|55|). Defining 

^ _ g2i drj ui 

= uo + ma, 



(59) 



eq. 



can be rewritten as 



/ + + Q^ + iQ)f = 0. 



(60) 



(61) 



An analogous equation has been derived in ref. [|T^]. Equation (|6TD is identical to the 
familiar equation for a spin-1/2 fermion in a time- varying background, if fl is identified 
with ma. To get a more transparent interpretation of this frequency Q, it is useful to 
consider the limit |$| ^ Mp, as suggested in ref. [0. Since m oc Mp^, in this limit the 
function G tends to the expression 

p — 2imMp 



G 



P 



(62) 



15 



where p = \d(^/dt\'^ + and p = \d^/dt\'^ - \dq,W\'^. One can now verify that G 

satisfies the following differential equation 

G d^W 

Therefore, in the limit of Mp —>■ oo and fixed |$|, u — »• —d^W and the frequency Q 
of the oscillations corresponds to the superpotential mass parameter of the Goldstino 



which is 'eaten' by the gravitino when supersymmetry is broken. Therefore, eq. (61) 
describing the production of helicity-1/2 gravitinos in supergravity reduces, in the limit 
of |$| <C Mp, to the equation describing the time evolution of the helicity-1/2 Goldstino 
in global supersymmetry and no suppression by powers of Mp is present. The frequency 
of the oscillations fl depends upon all the mass scales appearing in the problem, namely 
the Goldstino mass parameter d'^W/d^'^, the Hubble rate H and the gravitino mass m. 
The production of the helicity-1/2 gravitino is expected to be dominated by the fastest 
time-varying mass scale in the problem. 

To compute the abundance of helicity-1/2 gravitinos generated during and after 
an infiationary stage in the early Universe, one needs to discriminate among various 
supersymmetric infiationary models |l^. A crucial point to keep in mind is that a 
generic supersymmetric infiationary stage dominated by an F-term has the problem 
that the fiatness of the potential is spoiled by supergravity corrections or, in other 
words, the slow-roll parameter t] = MpV"/V gets contributions of order unity. In 
simple one chiral field models based on superpotentials of the type W = m$^/2 or 
higher powers in $, W ~ supergravity corrections make infiation impossible to 
start. To construct a model of infiation in the context of supergravity, one must either 



invoke accidental cancellations [|T^, or a period of infiation dominated by a D-teim ||18|| , 
or some particular properties based on string theory [|19|. 

Since the theory of production of helicity-1/2 gravitinos looks similar to the case 
of helicity-1/2 fermions with a frequency Q, one can use as a guide the recent results 
obtained in the theory of generation of Dirac fermions during preheating [^. During 
infiation, since the mass scales present in the frequency Q are approximately constant 
in time, one does not expect a significant production of gravitinos (the number density 
can be at most 113/2 ~ Hf, where Hj is the value of the Hubble rate during infiation). 
However, in the evolution of the Universe subsequent to infiation, a large amount of 
gravitinos can be produced. During the infiaton oscillations, the Fermi distribution 
function is rapidly saturated up to some maximum value of the momentum k, i.e — 
1 for k < /cmax and it is zero otherwise. The resulting number density is therefore 
^3/2 ~ ^max- The value of fcmax IS expectcd to be roughly of the order of the inverse of 
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the time-scale needed for the change of the mass scales of the problem at hand. Let us 
give a realistic example. Consider the superpotentialQ 

W = S{K^(j)- 12^), (64) 

where k is a dimensionless coupling of order unity [^|. Here, and are oppo- 
sitely charged under all symmetries so that their product is invariant. The canonically- 
normalized inflaton field is $ = -\/2|'S'|. The superpotential ( |64|) leads to hybrid inflation. 
Indeed, for $ ^ $c = f^/y/^, = = and the potential reduces to V = fi'^ plus su- 
pergravity and logarithmic corrections [17]. Therefore, in this regime the Universe is 
trapped in the false vacuum and we have slow-roll inflation. The scale fi is flxed to be 
around 10^^ GeV to reproduce the observed temperature anisotropy. Notice that in this 
period, the Goldstino mass d'^W/dS'^ is vanishing. 

When $ = inflation ends because the false vacuum becomes unstable and the 
fields and rapidly oscillate around the minimum of the potential at (00) = /i^/k, 
while the field $ rapidly oscillates around zero. The time-scale of the oscillations is 
0(/i^^). The mass scales at the end of infiation change by an amount of order of fi 
within a time-scale ~ fi"^. Therefore, one expects fcmax ~ and 71,3/2 ~ a/i^, where a 
summarizes the uncertainty in the estimate. If all the energy residing in the vacuum 
during infiation is istantaneously transferred to radiation, the reheating temperature 
would result to be Trh fi 10^^ GeV and the ratio 12^/2 to the entropy density would 
be ~ tt- This is not a realistic situation - though - because such a high reheat- 

ing temperature is already ruled out by considerations about the gravitino generation 
through thermal collisions, see the bound in eq. (^). In a more realistic scenario in which 
reheating and thermalization occur sufficiently late, the number density of gravitinos 
decreases as - a being the scale factor - in the post-infiationary scenario, presum- 
ably characterized by a matter-dominated Universe. If this is the case, at reheating the 
final ratio ^3/2 to the entropy density is 

^ ~ a^. (65) 
s fx 

If a = C(l), this violates the bound in eq. (|l|) by at least five orders of magnitude even 
if Trh ~ 10^ GeV and imposes a stringent upper bound on the reheating temperature 
Trh ~ 1 TeV. Notice that the non-thermal production is about five orders of magnitude 
more efficient than the generation through thermal scatterings during the reheating 
stage, irrespectively of the value of Trh- A similar result may be obtained in the case 
of D-term infiation. 



^Our results on gravitino production are based on a simple one chiral superfield model, but we 
expect similar results to be valid in the case of a theory with more than one superfield. 
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The ultimate reason for such a copious generation of gravitinos is that the system 
relaxes to the minimum in a time-scale much shorter than the Hubble time ~ H]'^, 
since the frequency is set by the height of the potential F^/^ » Hj during inflation. 
This is a common feature of realistic models of supersymmetric inflation. 

4 Conclusions 

In conclusion, we have investigated the non-thermal production in the early Universe 
of particles with mass in the TeV range and couplings to matter suppressed by powers 
of Mp. If produced with too large abundances, the late decays of these gravitational 
relics may jeopardize the successful predictions of standard big-bang nucleosynthesis. 
We have shown that scalar moduli - which are generically present in the mass spectrum 
of supergravity and string theories - may populate the Universe in large amounts as a 
result of the expansion of the background space-time acting on the vacuum quantum 
fluctuations of the moduli fields. The resulting number to entropy density ratio depends 
upon the way these scalar moduli are coupled to gravity and what is their effective mass 
during and after inflation. We have shown that the generation of these moduli fields 
poses no problem only if they couple to gravity with ^ close to 1/6, i.e. they are 
conformally coupled, and if the effective mass squared ChH"^ is comparable or larger 
than the Hubble rate, i.e Ch ~ 1- On the contrary, if ^ ^ 1/6 and Ch ~ 1, the 
generation can be so efficient that the standard predictions of nucleosynthesis are safe 
only if the final reheating temperature is as low as 100 GeV. This upper bound is 
considerably smaller than the bound of about 10^ GeV obtained from considerations 
about the production of scalar moduli through thermal collisions during the reheat 
stage. 

In this paper, we have also studied the non-thermal generation of gravitinos in the 
early Universe. The spin-3/2 part of the gravitino is excited only in very small amounts, 
as the resulting abundance is proportional to the small gravitino mass. On the other 
hand, the spin- 1/2 part obeys the equation of motion of a normal helicity-1/2 Dirac 
particle in a background whose frequency is a combination of the different mass scales 
at hand: the superpotential mass parameter of the fermionic superpartner of the scalar 
field whose F-term breaks supersymmetry, the Hubble rate and the gravitino mass. 
Again, in most realistic models of supersymmetric inflation, the non-thermal production 
of gravitinos turns out to be much more efficient than their thermal generation during 
the reheat stage after inflation. We have found, for example, that the ratio n3/2/s for 
hcIicity-1/2 gravitinos in generic models of inflation is roughly given by Trh/V^/'^., where 
~ 10^^ GeV is the height of the potential during inflation. This generically leads 



18 



to an overproduction of gravitinos. However, we should stress that the non-thermal 
gravitino production is quite sensitive to the specific infiationary model. Therefore, it 
will become an important ingredient in the search for realistic supersymmetric models 
of inflation. 
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